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Abstract 

A new method is presented for solving the momentum-space Schrodinger 

equation with a linear potential. The Lande-subtracted momentum space 

integral equation can be transformed into a matrix equation by the Nystrom 

method. The method produces only approximate eigenvalues in the cases of 

singular potentials such as the linear potential. The eigenvalues generated 

by the Nystrom method can be improved by calculating the numerical errors 

and adding the appropriate corrections. The end results are more accurate 

eigenvalues than those generated by the basis function method. The method is 

also shown to work for a relativistic equation such as the Thompson equation. 
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I. INTRODUCTION 



The momentum space Schrodinger equation has a singular kernel for both the Coulomb 
and linear potentials. The Coulomb singularity is removed with the Lande subtraction 
method JI],0]- Previous work [0-0] showed how to remove the singularity from the linear 
potential using a subtraction method with basis functions. A problem with this method is 
that one must guess a suitable set of basis functions in advance. In this paper, we show that 
Nystrom method J7| can solve the same problem more simply amd accurately. We begin 
with a review of the basis function method. Then we introduce the Nystrom method and 
apply it to the s-state momentum space Schrodinger equation with a linear potential. We 
use our new numerical results to show that the Nystrom plus correction method is more 
accurate than the basis function method. At the end, we generalize the Nystrom method to 
higher angular momentum quantum numbers (/ > 0). 



II. THE BASIS FUNCTION METHOD 

We begin this paper with a discussion of the basis function method to give the proper 
theoretical motivation. We shall use the simplest momentum space Schrodinger equation to 
illustrate the principles of the numerical methods, which is the s-state equation. 

The momentum space Schrodinger equation is related to an integral equation of the 
Fredholm type 

J K(p,p')(f)(p)dp = X(f)(p). (1) 
Suppose that the wavefunction can be expanded in a set of basis functions {(?«}, such that 

N 

<j>(p) = J2 C i9i(p), (2) 
i=i 

where Ci are constant coefficients. Substitute Eq. into Eq. |T|] to obtain 

N N 

£ / *(P. P') C i9i ip') dp' = \Y. CiQi (p) • (3) 

i=l J i=l 
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Now multiply both sides of Eq. || with gj(p) and integrate over p to symmetrize the equation 
over i and j, 

N „ „ N 



^2 C i[ J k (p,p') 9i{p')9j{p) dp' dp = A C i j 9i(p)9j(p) dp 
*=i "> . ' »=i "> .. " 



(4) 



and the result is a matrix equation, 



N N 

2_j A-ijCi = A BijCi 
i=i i=i 



(5) 



where C, is the eigenvector and A the eigenvalue. The indices i and j correspond to the 
quadrature points p and p' . N represents the number of mesh points. In the case of the 
momentum space Schrodinger equation with a Coulomb or linear potential, the kernel A 
is singular. A simple example is the momentum space Schrodinger equation with a linear 
potential in the S 1 - state pLH, 



f-</»no(p) + — \ [ 

2/J, TXJ) Jo 



A/ ' v 



■r-OS(y) + QM 
pp 



<Pno(p') dp' = E n0 (f) n0 (p), 



(6) 



v l (p,p') 



where y = (p 2 + p' 2 )/2p'p, 

QM = p'p 

and 

V 



(p + p') 2 + rf (p — p') 2 + rf 



(7) 



^Q'M = rf (p 2 + P' 2 + rf 
pp v 



(p + p') 2 + rf (p — p') 2 + T] 2 

Lande subtraction fHHH involves subtracting a zero term 



(8) 



^QM + Q'o(y) 
pp 



dp' = 



(9) 



from Eq. || such that 

— Mp) + — / 

Ifi irp z Jo 



pp 



[<t>no(p') ~ 4>no{p)] dp' = E n0 <p n o{p). (10) 

Using Eqs. the integral in Eq. JTD| for p > in the limit of y — > 1 can be shown to 

equal 



lim lim — 

rj— »0 p— »p' 7T 



1 



(P - P') 



A2 



dp 



0. 



(p — p') 2 + r] 2 J (p — p') 2 + r] 2 

The order of the limits in Eq. JTTJ] is important. The reverse order will lead to the nonsensical 
result / Q' (y) dp' = 0. Next, in the limit of p, p' — > 0, (p + p') 2 = (p — p) 2 . By substituting 
this equality into Eqs. it can be shown again that the integral in Eq. [TIJ vanishes for 
p — > at y = 1. At the end, the integral vanishes at y = 1, Vp. Away from the singularities, 
both integrands in the integral of Eq. Jl0| are finite. By taking 77 — > 0, the first integrand 
vanishes. The final form of Eq. JT(| is 

A; 



TtMP) + 2 , 

2/i 7r]r Jo 



Qd(l/) [0no(p') - 4>no{p)} dp' = E n0 (j) n o(p), 



(12) 



where Q' (y) = 1/(1 — y 2 ). As mentioned before, is expanded in basis functions, followed 
by integrating Eq. |12[ over p to generate a matrix equation. The basis functions used in 
previous publications are 



g f(p) = exp 



and 



gfip) 



—p 2 i 2 
M 



1 



(13) 



(14) 



(t/M) 2 +p 4 

respectively, where M is the maximum number of basis functions used. M has a maximum 
because the code crashes when too many basis functions are used. The basis functions gf(p) 
and gfip) have M = 18. 

The singularity of the kernel is a major challenge in solving the integral equation with 
a linear potential. It was shown 0] that a simple pole remains even after subtraction. The 
role of the basis functions is to make possible the evaluation of the Cauchy principal value 
of the subtracted integral using the Sloan method |J. To illustrate the Sloan method, we 
suppose that f(x) has a simple pole such that 

9{x) 



fix) 



x 



(15) 
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where g(x) is regular. The Cauchy principal value of the subtracted integral of f(x) can be 
evaluated if the range of integration is symmetric. For example, the numerical integration 
of 

/' 3{X) ~ 9(0) rfx (16) 

J-l X 

yields the Cauchy principal value because the point x = is skipped when quadrature points 
are generated in the symmetric interval (—1,1). The subtraction term has zero contribution 
since 

[ -dx=0. (17) 
7-i x 

The purpose of this term is to justify the existence of the Cauchy principle value and to 
reduce numerical errors. In order to apply the Sloan method to Eq. fl^fl , the integration 
variable is transformed from p to x such that x is centered at zero and its range is symmetric. 
In the case of the Coulomb potential, the kernel has a logarithmic singularity, 

(18) 

which is completely removed by Lande subtraction MM because no simple pole remains after 
the subtraction. 

The key to the success of the basis function method is the availability of a suitable set 
of basis functions for a particular problem. Unfortunately there is no a priori reason why 
the same set of basis functions will work in every situation. For this reason, it may be 
advantageous to have a method (such as the Nystrom method) that does not depend on a 
choice of the basis functions. 

III. THE NYSTROM METHOD 

In general, an integral equation of the Fredholm type 

G(pMp) + / F(p, p') 4>{p') dp' = A </>(p) 
Jo 
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Qo(y) = 2 ln 



(19) 



can be rewritten as a matrix equation as 



N N 

K V fa = H ( G H + fa = ^ fa' 

3=1 3=1 



(20) 



where Kij is the kernel and i and j are now indices corresponding to p and p' . Instead of 
integrating over p from to oo, we integrate over x from —1 to 1. Transform to Pi by the 
transformation 



p{.r) = tajj ( vr 



(21) 



The mesh points Xi and the weights wU are generated by the gaussian quadrature rule using 
the routine gauleg from Numerical Recipes @. In order to integrate along x G [—1,1] 
instead of p G [0, oo), Eq. |L9| is transformed as 



G(xU(x) + 1 1 F(x, x') 4>{x') dx' = X <j)(x) 
'-l ax' 



(22) 



Changing the dummy variable inside the integral and substituting the differentiation of 
Eq. with 



7T 2 ( \ + X \ , 7T / 

dp = — sec 
4 



I 7rJ drr = — (1 + p J dx 



(23) 



gives 



^fa+^J^Q'aiy^Vlfa-fa] ^(^pH dxi = E< 



(24) 



Eq. |24] can now be written matrix equation, 



P 



X N (\ + x- 

. 2/l ■ J^I J2Q'o(y ± 1 )<■>., sec 2 [—^ 7T ) ui, 



■j^fait,Qo(y + l ) se ° 2 ( ' 

^Pi k=l 
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7T ) wtu = E i 



(25) 



The left hand side of Eq. []25| is the kernel times the eigenvector and the right hand side 
is the eigenvalue times the eigenvector. The sum over k is independent of the eigenvector, 
which is just a scalar. The terms on the left that have only one index i belong to the diagonal 
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elements Ka. The terms with mixed indices % and j make up the off-diagonal elements, K^. 
More explicitly, the matrix elements of the kernel are 

K« = ~ ^ £ Q'oiv^ 1) (l + Pi) (26) 



t k 



K H = ^ Qo(v) (i + P$) (i * j)- (27) 

So far the kernel K is asymmetric under the interchange of % and j . We can improve the 
stability and the efficiency of the numerical solutions by symmetrizing Eq. |p5|| . We do so 



by multiplying the equation with p\ (1 It will change the original matrix equation 

Kx = Ax (28) 

to an equivalent matrix equation 

K'-x = AC-x, (29) 

where C is a diagonal matrix and K' = C ■ K. If C is positive definite, meaning 

x T • C • x > 0, V vectors x, (30) 

then C can be Cholesky-decomposed as 

C = L-L T , (31) 

where L is a unique lower triangular matrix. The reason for performing Cholesky decompo- 
sition is that the new matrix 

K" = L- 1 • K' ■ (L-y (32) 
is real symmetric and yields the same eigenvalues as Eq. |j2"9"| . In the case of Cu — pi (1 +pf), 



C is guaranteed to be positive definite. After symmetrization, Eq. [^TJ] does not change 



{K'li = Ka) while Eq. |27j becomes 



K "i = £^o(?/) JO- +Pf) (l +P]) ™tMv ii^j)- (33) 



The eigenvalues of K" can be calculated by using standard packages such as EISPAK. 
In this paper, we use the tred2 and tqli routines in Numerical Recipes f7|. 



IV. THE CORRECTION METHOD 



Maung, Norbury and Kahana have shown that the subtraction method does not 
completely remove the singularity at y = 1. There is a residual simple pole term 

f* ■ ■ if! (34) 

7T dp JO {p+P) [p—p) 

remaining after the subtraction. The basis function method evaluates the Cauchy principal 
value by the Sloan method as described in Section 2. The Sloan method eliminates the simple 
pole term by integrating symmetrically around the singularity. Symmetrical integration 
involves splitting the range of integration is into two intervals, 

dp' + 

tip 

The singularity at p — p' is contained in the first term on the right hand side of Eq. f35 



/•oo r2p roo 

/ dp' = dp' + / dp'. (35) 

JO JO J2p 



which is assigned a symmetric transformation rule (dp'/dx)\. The second term generally has 
a different transformation rule (dp 1 /dx)2 because it is mapping between two different sets, 
namely (2p, oo) — > (1, M] (for some real number M), such that 

Notice that the division of the range of integration moves with p. If two transformation 
rules are used with a moving division, each row (column) of the kernel has a different way 
to map [0, oo) to [—1, M]. But the eigenvector 0(p) must be mapped to fa in a unique way. 
This mismatch between the mappings of the kernel and the eigenvector does not affect the 
basis function method (see Eq. [2.24] of reference 0) 



N 



N 



/ 7T 9j(p)9i d P+ — / / Qo(y)9j(p) [9i(p') ~ 9i(p)] dp dp 

Jo III IT Jo Jo 



2(JL 

POO 

Ej^Ci P 2 9j(p)g i (p)dp (37) 

because the eigenvector C{ is an A^-tuple of the coefficients of the basis function expansion 
of the wavefunction <p(p) and is independent of the transformation rules. In the case of 



the Nystrom method, the problem is real, at least for the range of integration that we 
are interested in. Therefore we cannot evaluate the Cauchy principal value by symmetric 
integration in the Nystrom method. In other words, a new method is needed to treat the 
errors arising from the simple pole term. 

So far the error term Eq. J34J is not included in the Nystrom kernel in our derivation 



and is contributing to the errors of the eigenvalues. Since the error term Eq. |54j involves 
d(p/dp, we associate it with the error of the wavefunction 

where the mesh size Ap has an iV -1 dependence. This fact leads to an estimate of the N 
dependence of the error of the eigenvalue, AE. Let the approximate eigenvector to be <$' 
and the approximate eigenvalue E'. It is reasonable to say that an approximate kernel K 
acting on an approximate eigenvector <ft' yields an approximate eigenvalue E' as in 

K<f/ = E'(j)' (39) 
K((f) + A0) = (E + AE)((f) + A0). (40) 



It is easy to see that 



AE ~ (K - E) 
(K-E 



\ V J dp 

=4 < 41 > 

It is safe to assume that (K — E) d<p / 'dp « 1. can be interpreted as the normalization. 
The product of all of the pseudo-constants is labelled as the coefficient e. The approximate 
eigenvalue E' produced in the background of Eq. [R4| is related to the true eigenvalue E by 



E' n = E n + ef n (N), (42) 

where n is the principal quantum number, e a constant and f n {N) is a function approximately 
equals to iV -1 . In general, f n (N) varies slightly depending on the type of integral equation 
and the potential. As a first order approximation, assume that 

9 



f n (N) = N- 1 ~ a ^ 1 \ (43) 



The exponent of Eq. [^3| is a first order Taylor series expansion of some negative unity 



function around n — 1. The constant a is always taken to be small. More particularly, 
choose an a such that the variance of E n and e and \ 2 are minimized in the linear fit. 
Finally the refinement of an eigenvalue involves generating a set of E' n for various N by the 
Nystrom method and then extrapolating E n by a \ 2 linear fit in the graph of E n versus 
f n (N). In the case of Eq. a = 0.004 is an optimal choice. The numerical results are 



explained in Section 6. 

The order of the Nystrom algorithm is derived from those of tred2 and tqli, which is 
0(N 2 ) [0, compared with the basis function's 0(M 2 N), which comes from the product of 
the size of the matrix M 2 and the number of integration mesh points. N is typically around 
1000 and M is 20. The basis function method is generally more efficient than the Nystrom 
method. However, for any given set of basis functions, the accuracy of the eigenvalues cannot 
be improved arbitrarily by increasing the number of basis functions because M is bounded 
from above due to numerical errors. The prospect of improving the accuracy of the basis 
function algorithm depends on the availability of a set of more suitable basis functions for a 
specific problem. In the case of the Nystrom plus correction method, accuracy is optimized 
automatically by the correction scheme. The numerical results obtained by the Nystrom 
and basis function methods are quoted with optimal accuracy in this paper. 



V. EXACT S-STATE SOLUTION 



The eigenvalue of Eq. JT^| can be solved exactly in configuation space. We shall use the 
analytic results to check our numerical results. The non-relativistic Schrodinger equation 
can be written as 



^ + -4-) R-2^\ L r-E]R = 0. 
ar z r ar I 



(44) 



Let S = rR, then Eq. [@4| 

can be simplified as 
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Define a new variable 



-^S-2fj,[\ L r-E}S = 0. (45) 



x=(^Y[\ L r-E}, (46) 



such that Eq. |44| can be transformed as 

S" - xS = 0, (47) 

which is the Airy equation. The solution which satisfies the boundary condition S — > as 
x — > oo is the Airy function Ai(x). It is easy to show that the eigen-energy formula is 

where x n is the n-th zero of the Airy function counting from x = along — x. In reference [f|, 
the values Al = 5 and \i = 0.75 are used. In this case, the eigen-energy formula is 

E n = -2.554364772 x n . (49) 



VI. NUMERICAL RESULTS FOR THE S-STATE 

The accuracy of the Nystrom plus correction method is sensitive to the range of N. In 
this paper, increments of 100 in the range of 100 < N < 1400 are used. The reason for this 
choice is that there are not enough spacings between the eigenvalues for N < 100 and for 
N > 1500 the numerical noise begins to corrupt the monotonic convergent behavior of the 
eigenvalues. The correct eigenvalues are extrapolated from these numerical data by a x 2 
linear fit as described in Section 4. The exact S-state eigenvalues are tabulated against the 
numeical results obtained by the basis function method and the Nystrom plus corrections 
method in Table [||. It shows that the numerical results obtained by the Nystrom method 
plus corrections are more accurate than the results obtained by the basis function method. 
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The kernel written for the non-relativistic Schrodinger equation can be easily generalized 
to that of the relativistic 2-body Thompson equation in the center-of-mass frame by the 
replacement 

|- -> 2 (y^ 2 + m 2 - , (50) 

where p, is the reduced mass and m is the mass of each of the two equal mass elementary 
particles. The numerical results obtained using the Thompson equation is compared against 
those using the non-relativistic Schrodinger equation in Table § | calculating to 2 decimal 
places. Our new results are exactly the same as the previous results obtained in reference || 
that uses basis functions gf(p) from Eq. fT3H . 

VII. L^O KERNELS 

The / kernels for the linear and Coulomb potentials contain the Legendre function of 
the second kind Qi{y) and its derivative respectively. There are several mathematical issues 
that need to be addressed before constructing the I 7^ kernels. First of all, the definition 
of 

p 2 +p' 2 1 fp p'\ 
y =^p- = 2{^ + ^) (51) 

is easily seen to yield y > 1 for p,p r > 0. In reference ||, Maung et al. use the Legendre 
identity 



Qi{y) = Pi{y)Qo{y) - wi-i(y), 
1 1 

™i-i(y) = E -Pi-m(y)P m -i(y), (52) 



m=i m 



which is valid for —1 < y < 1 PJT0|] but can be extended to y > 1 by analytic contin- 
uation [ [Ll[| . Qi(y) is easily obtained by straightforward differentiation. The derivative of 
Legendre polynomial can be calculated from one of the recurrence formulas, 

-^- = y^y— + lp Uy), (53) 
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which can be computed numerically by a recursive call. The Legendre function can be 
generated by modifying the routine plgndr in Numerical Recipes J7| to allow y > 1. The 
accuracy of Eq. [|52[] and its derivatives are generally sufficient. Slightly more accurate results 
can be obtained by the explicit evaluation of the Neumann integral, 

G'frH/.'fv^W* (54) 



with derivative 



Q ' i{y) = ~\L^hf mdt - (55) 



The first few Qi(y) are 



Qod/) = ^ln^, (56) 
Q 1 (y) = iyln^±i-1, (57) 



Q 2 (2/) = i (3i/ i -lj ln^- I --y > (58) 

Q3(y) = i(5y 3 -3y)ln^±j-^ 2 + | (59) 

= ^ (35y 4 - 30y 2 + 3) In ^±L _ | y * + g y> (60 ) 

Q B (lO = ^ (63y 5 - 70y 3 + 15y) In j^±I - f y 4 + f y 2 - A. (61) 

<5K?/) can be obtained by the direct differentiation of Qi(y), such that 

Q'oiv) = y^- 2 , (62) 



y 

Q'M = -^— 2 - 1 -\n y —±, (03) 
1 - r 2 y + 1 

^^r^-^ 1 ' 1 ^!- 3 ' (64) 

Q > M = ^-^^ l Ay, ( 65 ) 

1 — y 4 y + 1 2 



„. , s 1 35w 3 — 15y , v — 1 35 9 5 . . 

^"T^—h-^^-T^ + i' (66) 



y 315y 4 - 210?/ 2 + 15 , y - 1 315 „ 105 

Q 5 (y) = In y —— - — y 3 + — y. (67) 

y 16 y + 1 
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As y — > oo, it is easily seen that Qo(y) = Qb(y) — > 0. This limit is true for all Qi(y) and 
Q'i(y) from applying the L'Hopital rule. Unfortunately straightforward numerical calculation 
of Qi(y) and Q\(y) by using Eqs. |56| - |67| leads to serious numerical errors as y — > oo. At 
the same time, it is observed that the numerical integration of Eq. |54| , |55|| are reasonably 
accurate in the same regime. Therefore the two representations are combined to minimize 
numerical error by using the Neumann integrals for y > y Q and the explicit formulas for 
V < Vo- Our codes use yo = 50. 

The subtracted momentum space NRSE with a linear potential is given in reference [[|, 
which can be simplified as 



^Mp) + — 2 T QKy) Mp) dp' - — 2 Mp) I QM dp' 

Zfl irp^ Jo 7rp Jo 

\ L l(l+l) f ooQ \ L7r 1(1 + 1) 



<Pm (p) r ^dp> + ^ <p nl (p) = E nl <p nl (p) . (68) 

Jo V V 4 



■up I Jo p p 

The matrix elements of a symmetric kernel for arbitrary / are 

= 73 EGofo^i) (i+ri) 



-^t^(l) (l+ft 2 ) (69) 

= j^- Qi(y) JO- +p") (i +pj) ^i^v i^j)- (70) 

QPiPj " 

Despite our method to control numerical noise, numerical errors still manifest themselves in 
the form of spurious large negative eigenvalues for I > 8. Fortunately the rest of the positive 
eigenvalues are accurate. Some sample eigenvalues for < / < 5 are shown in Table |TI 



which also compares the eigenvalues generated by both the p-space and r-space codes. The 
r-space eigenvalues are calculated by solving NRSE using the relaxation method |T3| . 

The Lande-subtraced momentum space NRSE equation with a Coulomb potential is also 
given in reference 0] and is simplified as 

ttMp) + — / Qi(y) <t>rdip') P dp' 

III Tip Jo 
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-p<j>rd(p) I -^T~dp' + -^p<j) n l(p) = E n 4nl{p)- (71) 

71 Jo p 2 

The kernel of a Coulomb potential can be symmetrized the same way as that of a linear 

potential. The matrix elements are 

V Pi X C \ - Qo{y + 1) (, . 2 \ , , ^CXPi 

(i+pi) Mi (72) 
Kij = ^fQi(y))/0-+ Pi) (l + Pfj wti wtj, {1+3). (73) 

The correction method which we have developed for the linear potential cannot be used in 
the Coulomb case. The only available technique of refining the eigenvalues of a Coulomb 
potential is by the way of increasing the number of mesh steps N. Some sample eigenvalues 
are shown in Table |IV|]. Since both the linear and Coulomb potentials can be symmetrized 
using the same formalism, we can easily splice the two kernels together to calculate the 
eigenvalues of the Cornell (linear plus Coulomb) potential 

V(r) = ^ + X L r. (74) 

It is not surprising that the correction method derived for the linear potential may also 
work for the Cornell potential because we expect that the error of the Cornell potential is 
dominated by the error of the linear potential term. But it is a surprise that the correc- 
tion method works more accurately with the Cornell potential than the linear potential as 
evidenced by vanishingly small variance and x 2 - 

VIII. CONCLUSION 

The basis function method requires a priori knowledge of the eigenf unctions in order to 
pick out an appropriate set of basis functions. The advantage of the Nystrom method is 
that no such prior knowledge of the eigenfunctions is needed. The kernel constructed by the 
Nystrom method is also much simpler than that by the basis function method. The eigen- 
functions can be generated by the same Nystrom routines that compute the eigenvalues. 
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The Nystrom plus correction is more accurate than the basis function method in the cases 
studied in this paper. In other words, the new method has all of the advantages-elegance, 
accuracy and versatility. In addition, the kernel of the relativisitc and non-relativistic equa- 
tion of motion with the Coulomb and linear potential can be symmetrized in exactly the 
same manner. It allows the calculation of the eigenvalues of a Cornell potential readily. Since 
the Nystrom method can be generalized for higher Vs, we can use it to calculate the Regge 
trajectories. Since the Thompson equation which we have solved is a 2-body equation, we 



can use it to analyze the experimental meson Regge trajectories [|1J]. This will be pursued 
in later work. 
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TABLES 

TABLE I. Comparisons of eigen-energies in GeV of the non-relativistic Schrodinger equation 
with a linear potential between the Nystrom method and the basis Function (BF) method. The 
basis functions being referred to here are gf(p) = \{i/M) 2 +p A ]~ l . The values of I = 0, X L = 5 GeV 
and ji = 0.75 GeV are used. 



n 


N = 100 


Nystrom 
N = 700 


N = 1400 


Corrected 


BF 


Exact 


1 


5.899211 


5.961921 


5.967339 


5.972379 


5.972 


5.972379 


2 


10.268443 


10.417386 


10.430047 


10.442010 


10.443 


10.442114 


3 


13.767781 


14.054263 


14.078517 


14.101276 


14.104 


14.101524 


4 


16.784747 


17.258395 


17.297500 


17.335360 


17.335 


17.335728 


5 


19.467512 


20.177458 


20.234722 


20.291708 


20.293 


20.292215 


6 


21.891635 


22.887999 


22.967933 


23.046820 


23.053 


23.047142 


7 


24.101339 


25.435892 


25.541743 


25.646532 


25.648 


25.646268 


8 


26.124257 


27.851711 


27.986463 


28.121481 


27.947 


28.120787 


9 


27.977844 


30.156480 


30.323418 


30.493311 


30.194 


30.488938 


10 


29.672260 


32.366010 


32.568895 


32.778297 


33.340 


32.769375 
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TABLE II. Comparisons of the ratios of Eigen-energies E n+ \/E\ using the Thompson equation 
(TE) and the non-relativistic Schrodinger equation (NRSE) using I = and A L = 0.2 GeV 2 . Mass 



is measured in GeV. 

n TE NRSE Mass 

1 1.72 1.75 1.5 

2 2.30 2.36 1.5 

3 2.80 2.90 1.5 

1 1.67 1.75 0.5 

2 2.18 2.36 0.5 

3 2.62 2.90 0.5 

1 1.63 1.75 0.3 

2 2.11 2.36 0.3 

3 2.51 2.90 0.3 
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TABLE III. Eigen-energies in GeV of the non-relativistic Schrodinger equation in momen- 
tum space (pNRSE) compared with those in the configuration space (rNRSE) and the relativistic 
Thompson equation in momentum space (TE). The r-space Thompson equation result is not avail- 
able. The values of n = 1, Xl = 5 GeV and n = 0.75 GeV are used. 







piN i v.j 






i in i v.i j 


7 

t 


at i nn 
iv = 1UU 


iv = /UU 


at i/i nn 
iv = 14UU 


Corrected 


Approx. 


n 
u 


^ 8QQ91 1 


^ Qfi1 Q91 

o.yoiyzr 


o.yo i ooy 


^ 079^70. 

o.y i zo i y 


^ Q71 Q 

o.y ( iy 


i 
l 


O.OZO ( zo 


o.o i/ilo 


o.0oz41o 


o.oobUUz 


o.oooU 


o 


1 D 89^nQQ 


1 D 847^3 


1 D 84Qfi7 c > 


1 D 8^1 'Wfi 


i n 8^1 4 


Q 
O 


1 9 Q1 71 94 


1 9 QD4991 


1 9 QD981 ^ 


1 9 QD91 1 7 
±z.yuz±± i 


1 9 QD90, 
iz.yuzu 


A 


14 874948 


14 81 9499 


14 8 nc i4fi9 


14 8D1 ^8 


1 4 Q7QD 
i-^i. y ( yu 


J 


1 7^n^&^ 

1U. 1 ouoou 


lu.UUUUOl 


iu.oy / uou 


1 ft ^Sfi^fil 

lU.dOUOUl 


1 ft ^84^ 






TT? 

J. Hj 








I 


V = 100 


iV = 700 


TV = 1400 


Corrected 


NA 





5.859885 


5.914287 


5.919054 


5.923117 




1 


8.164379 


8.202282 


8.205185 


8.208610 




2 


10.053574 


10.067261 


10.068464 


10.069762 




3 


11.700322 


11.680163 


11.678063 


11.676817 




4 


13.185124 


13.121767 


13.116634 


13.111239 




5 


14.553134 


14.437612 


14.427702 


14.418173 
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TABLE IV. Eigen-energies in eV of the hydrogen atom according to the non-relativistic 
Schrodinger equation with n = 1. 

I N = 100 N = 1400 N = 3000 Exact 

-25.286631 -13.600349 -13.598508 -13.598289 

1 -4.579043 -3.400415 -3.399659 -3.399572 

2 -1.463504 -1.511499 -1.510980 -1.510921 

3 -0.634523 -0.850358 -0.849940 -0.849893 

4 -0.329730 -0.544332 -0.543972 -0.543932 
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